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Using a metric perturbation method, we study gravitational waves from a test particle scattered 
by a spherically symmetric relativistic star. We calculate the energy spectrum and the waveform 
of gravitational waves for axial modes. Since metric perturbations in axial modes do not couple to 
the matter fluid of the star, emitted waves for a normal neutron star show only one peak in the 
spectrum, which corresponds to the orbital frequency at the turning point, where the gravitational 
field is strongest. However, for an ultracompact star (the radius R < 3M), another type of resonant 
' periodic peak appears in the spectrum. This is just because of an excitation by a scattered particle 

of axial quasinormal modes, which were found by Chandrasekhar and Ferrari. This excitation comes 
from the existence of the potential minimum inside of a star. We also find for an ultracompact star 
many small periodic peaks at the frequency region beyond the maximum of the potential, which 
^ , would be due to a resonance of two waves reflected by two potential barriers (Regge- Wheeler type 

and one at the center of the star). Such resonant peaks appear neither for a normal neutron star nor 
for a Schwarzschild black hole. Consequently, even if we analyze the energy spectrum of gravitational 
waves only for axial modes, it would be possible to distinguish between an ultracompact star and a 
normal neutron star (or a Schwarzschild black hole). 
PACS number(s): 04.25.Nx, 04.30.-w, 04.40. Dg 
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A laser interferometer to detect gravitational waves, such as the Laser Interferometric Gravitational Wave Obser- 
vatory(LIGO), VIRGO, GEO600, or TAMA300, will be in operation within several years [Q. It is one of the most 
urgent subjects for theoretical relativists to make a set of templates of gravitational waves. The direct detection of 
gravitational waves is very important not only as a new "eye" observing the Universe (gravitational wave astronomy) 



O 

but also as a new probe for fundamental physics. For example, we know that a neutron star, which is observed as a 
^Jy pulsar, has a mass M » IAMq and a radius R w 10 km. Although we can guess its central region with some theoret- 
ical ansatz, we have no direct information from the inside of a neutron star because neither radiation nor a neutrino 
is transparent in the central region. The direct detection of gravitational waves may provide some information about 
the inside of a neutron star when it is formed, resulting in new observational facts or constraints about the equation 
■ of state for high density matter. 

One of the most promising astronomical sources for these ground-based detectors is a coalesced compact binary 
such as a neutron star binary. To complete a set of templates from such a source, many studies about the emitted 
gravitational waves have been done using various techniques. 

A numerical simulation without approximations (numerical relativity) is probably one of the best ways to calculate 
gravitational waves emitted from the compact binary system |^,^)- Although steady progress in numerical relativity 
has been achieved, there are still some difficulties which prevent us from finding the final results. 

Instead, there are some powerful approximation methods to mimic the coalescence of a compact binary system. 
One such approximation is the perturbative approach. Regge and Wheeler Q and Zerilli first formulated a method 
for metric perturbations in a Schwarzschild black hole spacetime. Then there are many works using such a black hole 
perturbation method in which we describe gravitational waves by metric perturbations in a black hole background 
spacetime and treat a companion of the binary as a test particle |^ |8|. This approximation for gravitational waves 
in a head-on collision of two black holes gives a good agreement with the results obtained by numerical simulation 
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P). Thus, such a perturbative approach may provide a good approximation for gravitational waves from a compact 
binary as well. 

However, there are so far few works for neutron star perturbations. One may think that if we have two 1.4Mq 
neutron stars, it will be a black hole after coalescence, when the ringing tail of the gravitational waves can be described 
by quasinormal modes of the black hole. The quasinormal modes of a neutron star may have nothing to do with the 
final stage of the coalescence. However, before black hole formation, we have to analyze the orbital evolution of the 
neutron star binary. Therefore, we believe that a perturbation analysis for a neutron star is also important. 

For a nonradial pulsation of a spherically symmetric relativistic star, perturbation equations both for axial and 
polar modes were first derived by Thorne and Campolattaro |uj. Similar to the black hole perturbation, the axial 
mode for a spherical star is described by a second order differential equation. On the other hand, the polar mode was 
first described by a fifth order differential equation [^0|. However, because of the dynamical degrees of freedom of 
polar modes, i.e., sound waves and gravitational waves, the basic equation can be reduced to a fourth order differential 
equation as shown by Lindblom and Detweiler JjTfl . Lindblom, Mendell, and Ipser [jl2| also derived a couple of second 
order differential equations, which could be interpreted as those for sound waves and gravitational waves, respectively. 

As in the case of black hole perturbations, quasinormal modes, if they exist, will be important in the analysis 
of gravitational waves from a compact binary. For the polar mode, there are two types of quasinormal modes of a 
neutron star; one is a fluid oscillation mode similar to the /, p , g modes of a Newtonian star (cf. JT^j); the other is 
a wave mode (u> mode), which exists only in general relativity pTHlq]. Even for the axial mode, Chandrasekhar and 
Ferrari illustrated that quasinormal modes exist if the effective potential has a minimal as that in an ultracompact 
star. Kokkotas |20| numerically calculated such axial quasinormal modes for various ultracompact stars. Kojima |2l]] 
extended it to a slowly rotating star, setting Regge- Wheeler gauge conditions. 

As for emitted gravitational waves from a binary system, Kojima |22| studied them for the case that a test particle 
is moving in a circular orbit around a polytropic star. For a circular orbit, the axial mode cancels from the orbital 
symmetry so that the polar mode was only taken into account. He showed that the resonance mode appears in the 
energy flux when the orbit reaches the radius where the orbital frequency coincides with that of the quasinormal 
mode. For the axial mode, recently, Borrelli |23) calculated energy spectra of gravitational waves emitted by a test 
particle spiraling into an ultracompact star. In his result, its energy spectrum shows many peaks, which correspond 
to axial quasinormal modes of the star. Although it is important, since we are interested in a compact binary, a direct 
collision might show some different features from the binary case. We also do not know what will happen when a 
test particle reaches the surface of the star. To avoid such an unknown factor and to analyze the case of a compact 
binary, we first study gravitational waves emitted when a test particle is scattered by a neutron star in this paper. 

For a Schwarzschild black hole, Oohara and Nakamura § analyzed it and showed that the energy spectrum has 
no peak except for one which corresponds to the orbital frequency. Here we consider a spherical star instead of a 
Schwarzschild black hole. We study only the axial mode which does not exist in Newton gravity. To discuss the 
dependence of the equations of state, we analyze two models: a uniform density star as in jl^,|2^,^3j and a polytropic 
star. 

This paper is organized as follows. In Sees. [n|A and^B, we briefly review the perturbation theory for a spherically 
symmetric relativistic star. We show our numerical results in Sec. [II. Section ^ is devoted to a discussion. Some 



numerical techniques for calculating gravitational waves are summarized in the Appendix. Throughout this paper, 
we use units of c = G = 1 and a metric signature of (— , +, +, +). 

II. LINEARIZED EINSTEIN EQUATIONS: AXIAL MODES 



We start with summarizing briefly the metric perturbation of a spherically symmetric relativistic star. A spherically 
symmetric background metric gffi is described as 

ds 2 = gWdaPdaf = -e^dt 2 + e x ^dr 2 + r 2 (d6 2 + sin 2 6W0 2 ) , (2.1) 



e-*> = l-^, (2 . 2 ) 
r 

where M(r) is a mass function inside a radius r. We assume a perfect fluid as the stellar matter, i.e. 

T» v = {p + P)u^u u +Pg ttv , (2.3) 

where p and P are the total energy density and the pressure, respectively, which satisfy a barotropic equation of 
state, P = P(p). Using Tolman-Oppcnhcimer-Volkov equations, we construct a stellar model and its spacetime 
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metric numerically except for the constant density star case. The radius R of a star is defined by a surface where the 
pressure vanishes, and the mass of the star M is given by M = M(R). 

For linearized Einstein equations, we can decompose the perturbed metric into an axial mode h^* lal ^ and a polar 
mode /i^ olar ^: 

9„u = 9$ + + h% o1 ^. (2.4) 

Using Regge- Wheeler gauge conditions, we will integrate the interior and exterior regions of the linearized Einstein 
equations, separately, and impose appropriate junction conditions at the surface R to find the solution. In this paper, 
we only deal with axial mode perturbation. Then, we write down the linearized Einstein equation only for axial modes 
in the following subsections. 



A. Interior region of the star 

In order to find the linearized Einstein equations for axial modes, we expand the perturbed metric by tensor 
harmonics for polar angle 6 and azimuthal angle tfi > and proceed to a Fourier transformation for the time coordinate 
t. Although there are ten independent components for the linearized metric, four of them are gauge freedoms of 
coordinate transformations, which are fixed by a choice of Regge- Wheeler gauge 0, resulting in two components 
belonging to axial modes, while the rest of the four correspond to polar modes. For axial mode perturbations, we can 
set 

h { ^ iaX) dx' l dx v = -!- / dwV2/i 0! (r) (- , 1 - d ^ lm dtdO + sin O^^dtdA e~ wt 
F 2tt J ' \ smf) dm oO I 

l,m 
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where hoj and hu are functions of the radial coordinate r defined below, and Yi m (9,4>) is a spherical harmonics. 

The axial mode in the interior region of a star is described by a single wave equation for gravitational waves, because 
a rotation does not exist in a spherically symmetric background spacetime. That is, we have a single wave function 
Xj™^ , which determines the behaviors of axial modes as 

h 1>lu = e -(^)/ 2 rX£ nt) , (2.6) 

V*. = -J V - X) ' 2 i Kf t} ) • (2-7) 
The equation for the wave function x\™^ is derived from the linearized Einstein equation for axial modes as 



d 2 X^ t] 
dr* 2 



+ (uj 2 - y/ int) ) x£ nt) = 0, (2.8) 



where the "tortoise" coordinate r* and the effective potential V^ 1 "*"* are defined as 



e 



,(int)_ u ( + 6Af(r) 



Vi ' = ^ ( ~^~~2 ~ — -P ~ MP-P) ) > CI. 1.0 ) 

respectively. 

As for the 
around the center r = 0, we find 



As for the boundary condition, the wave function A ; ^ nt ^ must be regular at the center of a star. Expanding X^ 



-( int ) — ^. J J+i _i L 



2(21 + 3) 



4m/ + 2i ( i(2Z-l)p c -P c ) "■ 



r l+3 + 0(r l+5 ) } , (2.11) 



where r]^ is an arbitrary constant, and p c , P c , and v c are the central density, pressure, and metric functions, respec- 
tively. 

Apart from a free parameter rji^, we obtain a solution A^ 11 ''' by integrating Eq. (2.S) with the boundary condition 
( |2.11 ). Fu rther details of the numerical techniques constructing the wave function xj^ are given in the Appendix, 
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B. Exterior region of the star 



The background metric outside the star is the Schwarzschild spacetime, i.e. , Eq. (2.1), with 

2M 



e u = e-* = l . (2.12) 

r 

As for the trajectory of a test particle, without loss of generality, we can assume that the particle moves on the 
equatorial plane (9 = tt/2). Then, the equation of motion for the particle is described as 

(2.13) 

dR\~ ~ 9 / 2M\ I L 2 \ . . 

3F =E VS) 1 + * ■ (214) 
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(2.15) 
(2.16) 

where the particle trajectory is = (T(r), -R(t), O(r), $(t)), and ^t, _E, and L are the mass, the normalized energy, 
and the normalized orbital angular momentum of the particle, respectively, with E = E//i and L = L/fi. The 
normalized effective potential energy tAP artlclc ) (r) of the particle is given as 



y (particle) 



2M\ ( I 2 



r 2 



(2.17) 



In this paper, since we only consider the motion of a test particle scattered by a spherical star, the particle orbit is 
unbounded, i.e., 



L\/M>A, 1 < E < Vi P a T iclc \ (2-18) 
i of the potential. 

The perturbed metric outside a star only comes from an oscillation of spacetime, that is, gravitational waves. Then, 



where V^ax rtlcl °' ) is the maximum value of the potential 



similar to the interior region, the linearized Einstein equations for the metric /ii,; m and ho t i m defined as Eq. (2.5) are 



reduced to a single wave equation for a new wave function Xj^ defined by 

hi,imu=e x rX^\ (2.19) 
i -\ d ( „i cxt )\ 87ri r 2 e~ x 



rX icxt) 



h .i m u = [ rX\^J ) - — 7 ====D lm (u, r), (2.20) 



where n = (I — + 2)/2 and Di m (u>, r) is one of the source terms given below. 



From the linearized vacuum Einstein equations, we find the wave equation for X^^ (Regge- Wheeler equation) 



as 



, 2 v(ext) 



dr* 2 

where r* is the tortoise coordinate defined by 



r* = r + 2Mln(-^-l), (2.22) 



2M 



V^ oxt ' is Regge- Wheeler potential given by 



if->-.-»(!2£i-5£), ( , 23) 
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and <5;ff* is the source term, which is described by the energy-momentum tensor of a test particle expanded by tensor 



harmonics |j] as 



-2A 



v/n+1 



Qim(w, r) + 



^2n(n + 1) dr 



with Qi m (t,r) and D tm (t,r) being 

Qim(t,r) -- 



1 



Di m (t,r) 



y/n+1 r 

J_dYlrnd& 

sin6 <9<l dr 
1 



dT^-6(t ~ f(T))5(r ~ R(t)) 
dr 



^J2n(n+ 1) 



— sin 6 — ^ r- 

dT«(t - f (r))«(r - ii(r)) 



de 



sin 2 9 



dr" 



1 



sin 



s dO d$ - /As 



(2.24) 



(2.25) 



(2.26) 



Xirn, Wi m are the tensorial part of a tensor harmonics defined by 

d ( d 



Xl m — 2- 



m 



cot6> F im , 



1 d 2 



(2.27) 
(2.28) 



and the overbar denotes the complex conjugate. 

The boundary condition outside the star should be imposed such that there is no incoming wave at infinity, 



y(ext) A j, 



(r* -> oo), 



(2.29) 



where A; mu is the amplitude of an outgoing wave at infinity. The details of constructing the wave function X} are 
given in the Appendix, subsections A 2 and A 3 . 



C. Matching of two wave functions and gravitational waves 

We have to combine two wave functions, the interior wave function x/ mt ' ) and the exterior one X^ t ), at the surface 
of the star. The matching condition is that the wave function must be continuous and smooth at the surface. For the 
case of the axial mode, since the wave function inside the star does not couple to perturbations of the matter fluid, 
the condition turns out to be very simple as 

X^\R*) = X^(R*), (2.30) 
' (R*)=^^(R*), (2.31) 



dr* dr 
where R* is the radius of the star by use of the "tortoise" coordinate. 



Then we calculate the amplitude of th e gra vitati onal w aves Ai muJ as follows: We first integrate Eqs. ( |2.S| ) and (2.21) 



using their boundary conditions, Eqs. ( [2.11 ) and (2.2E). Next, we connect those solutions at the surface of the star 
using Eqs. (2.30) and (2.31), which determine two unknown coefficients f]i u and Ai mu} . 

From the wave function obtained as above, we obtain some information about gravitational waves emitted by a 
scattered particle. When we discuss the gravitational waves, we usually decompose them into two modes: + and x 
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modes. Then, we have to transform our coordinate system from the Regge- Wheeler gauge to the transverse-traceless 
one, which metric near infinity is described as 

ds 2 = ~dt 2 + dr 2 + r 2 (l + h + )d9 2 + r 2 sin 2 6(1 - h+)d<t> 2 + r 2 sia9h x d9d(j>. (2.32) 

Then the axial mode in transverse-traceless gauge at infinity turns out to be 



1 r°° rl 00 

±ih x =T— / —V V A lnu Vn(n + l)e iu ^-^ ±2 Y lm (d, <f>), (2.33) 



~°° " 1=1 m=-l 

where the spin-weighted spherical harmonics ±iXi m is defined by 



±2Yi m (6, <j>) = . 1 [W lm ± -4-X im ) . (2.34) 



i ; r- I ' ' flit -•— . p, 

2^Jn(n + 1) V sm& 
The energy-momentum tensor of the gravitational waves at infinity is given by 

T ^ W) = ^ (Re(V^S T) )Re(V^^( TT ))) , (2.35) 



32tt \ 

in the transverse-traceless gauge, where the angular brackets denote an appropriate average. 

From this definition, we find the energy spectrum dEg^ /du> of gravitational waves at infinity for the axial mode as 

^ff^^E E l(l + m-W + 2)\Ai™\ 2 . (2.36) 

1=2 rn=-l 



III. GRAVITATIONAL WAVES FROM A SCATTERED TEST PARTICLE 

A. Uniform density star 

In this subsection, we consider a uniform density star. We restrict our analysis only to the case of the / = 2 mode, 
just for simplicity. We find that features of the energy spectrum of emitted gravitational waves can be classified into 
two cases, which mainly depend on the shape of the effective potential. 

Figure [l] shows the effective potential [Eqs. (2.1C ) and ( J2.23 ) ] for a uniform density star with three different 



compactness parameters R/M. The effective potential has the following three features. First, the Regge- Wheeler 



potential [Eq. (2.23)] has a maximum point near r ~ 3M (accurately at r = 3.28Af for I = 2). Second, the potential 
is discontinuous at the surface because of a discontinuity of the energy density of the star at the surface. Finally, a 
potential minimum appears inside the star, if the star is ultracompact (for the case of R < 3M). 

According to the last feature, we study two cases separately. The first stellar type is the case that the surface is 
outside the potential barrier (R > 3M). For a usual neutron star, R ~ 5M, and then the first type is most likely 
a realistic astronomical object. The second stellar type is the case that the surface is inside the potential barrier 
(R < 3M), which may be an exotic star, if it exists. 

First we show our result for the first type. For the case of R — 5.0M, we show the energy spectra of emitted 
gravitational waves in Fig. ||. In this figure, there exists only a single peak in each spectrum. This peak corresponds 
to the orbital frequency of a test particle at the turning point, i.e., periastron. This means that most gravitational 
waves are emitted around the turning point, where the effect of the gravitational field is the strongest. No other feature 
is found in this spectrum. In fact, fixing the energy E and orbital angular momentum L, the spectrum depends less 
on the compactness parameter R/M (see Fig. ||). For a black hole background, we find quite similar behaviors in the 
energy spectrum (Fig. |^), except for the high frequency region, where the difference of the boundary conditions is 
crucial. We also find that the spectrum approaches closer to that of a black hole as the star gets smaller. 

We also calculate the waveform for the case of R = 5.0M and black holes, respectively. We see that a burst wave 
due to the encounter of a test particle emerges for each case. As the same as the energy spectrum, the difference 
between these waveform is very little (see Fig. ||) 

Therefore, for the first stellar type of R > 3M, the energy spectrum largely depends on the trajectory, but not on 
the background object. 
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Next we discuss an ultracompact star, which shows much more various structure in the spectrum. In this case, a 
minimum point of the effective potential appears inside the star. As for the trajectory of a test particle, we have two 
possibilities: one is that the particle gets into the inside region beyond the potential barrier and the other is that one 
passes through outside the barrier. Since the particle approaches closer to the star, it feels a stronger gravitational 
field. As a result, the excitation by the particle is more efficient and the emitted gravitational waves are enhanced. 
As a numerical example, we show those energy spectra for the model with radius R = 2.26M with the turning points 
# min = 3.145M and 6.628M in Fig. § 

Apart from a global maximum, which appears in the first type model as well and corresponds to the orbital 
frequency of a test particle at the turning point, we also find many small periodic peaks in Fig. ra. This microscopic 



peak structure is classified into two types, which are divided at the frequency of the maximum w max — yV l 



(~ 0.389 A/ -1 for I = 2) of the effective potential. In the lower frequency region (u> < ui max ), small and sharp peaks 
correspond to the quasinormal modes, as we will show it next. For the higher frequency region (u> > w max ), we find 
small periodic peaks, which may appear because of a resonance between two potential barriers. We will also discuss 
this later. 

In order to analyze our numerical results, we first have to know the axial quasinormal modes. Our method to 
calculate a quasinormal mode is based on the continued fraction expansion method which was originally used for 
black holes by Leaver |24j], and adopted for the polar mode of a spherical star by Leins, Nollcrt, and Soffel p6[ . We 
apply this method to axial modes as well. We show our results in Table | with those of Kokkotas [^(| whose method 
is different from ours. The table shows that both results agree quite well. 

In Fig. |^, we present both the quasinormal mode and energy spectrum. As seen from Fig. ^, axial "quasinormal" 
modes are classified into two: the lower frequency mode (to < w max ), which imaginary part is very small lm(uj) < 
O(I0~ 8 ), and the higher frequency mode (uj > w max ), which imaginary part is rather large (10 -8 ) < Im(tj) < 0(1). 
The reason is very simple. For a frequency lower than w max , we can have a quasibound state. Since trapped waves 
tunnel through the potential barrier, the energy will decrease. It is similar to the quasinormal mode of a black hole. 
Then, we may call it the axial quasinormal mode as well. However, if the frequency is higher than w max , we have 
no bound state. Nevertheless, we can find the "quasinormal" mode, by which we mean a wave solution with the 
same boundary condition as that of the conventional quasinormal mode, which is the outgoing wave condition. In the 
present case, although we find such wave solutions, those waves decay very quickly because the frequency is higher 
than the potential barrier. This is the reason why we find a large imaginary part of the mode frequency. The solutions 
which satisfy the boundary condition appear periodically as seen from Fig. |^. 

We find that the frequencies of these quasinormal modes coincide with those of spectrum peaks in the frequency 
region lower than w max . This is because a test particle excites those quasinormal modes. In fact, when we see the 
waveform (Fig. 0), the wave consists of two parts: large one, which corresponds to the burst wave by the encounter 
by a test particle, and the other part with several oscillations, which comes from exited quasinormal modes in the 
spectrum. Actually the waveform in the region after t — r* ~ 100M in Fig. [7] is similar to Kokkotas's result for a 
Gaussian pulse wave (2^] . 

On the other hand, in a frequency region higher than o; max , although many periodic peaks appear, those positions 
are not exactly the same as the frequencies of "quasinormal" modes. We may understand this fact as follows. Even 
if the wave frequency oj is higher than the potential barrier, the reflection coefficient does not vanish and oscillates 
with respect to oj, which period is determined by the potential form. In the present case, we have also an infinite 
wall near the center of a star. We then expect that small periodic peaks appearing in the spectrum are due to an 
interference between the incoming and outgoing waves reflected by two potential barriers. Since these waves decay 
very quickly as mentioned above, the process is quite dynamical, and then the peaks in the spectrum do not coincide 
with the frequencies of "quasinormal" modes. 

In Fig. ||, we show the energy spectra for two stellar models with R = 2.26M, and 2.3M, which correspond to the 
cases that the turning point is inside and outside the barrier, respectively. For fixed energy E and angular momentum 
L, the global peaks corresponding to the orbital frequency at the turning point agree with each other. On the other 
hand, the other small periodic peaks from the excited quasinormal modes do not agree with and rather depend on 
the compactness parameter R/M. Consequently, from the analysis of the energy spectrum of axial modes, we may 
determine the compactness parameter R/M of an ultracompact star. 

We can also compare a black hole and an ultracompact star. In Fig. ||, we show the energy spectra both for a black 
hole Q and an ultracompact star with R = 2.26M. From Fig. we see that even for the case of a Schwarzschild 
black hole, the global peak due to the orbital motion of the particle coincides with that for the star. However, unlike 
the stellar model, no small periodic peaks appear in the case of a black hole. In fact for a black hole, there is no 
peak after t — r*^ 100M (see Fig. 7). This is because the effective potential of the star has an infinite potential wall 
near the center, which guarantees the existence of quasinormal modes, while that of a Schwarzschild black hole has 
no such wall, and incoming waves fall entirely into the event horizon (ru = 2M). Hence, we can also distinguish a 
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background object, i.e., whether it is an ultracompact star or a black hole, analyzing only axial modes. 
As for the total energy of emitted gravitational waves £qw> we find 



e gw = e jrr, (3.1) 

with e = 2.501 x 1CT 2 and 3.541 x 1(T 2 for R = 5M and AM, respectively, and for E = 2/z and L = 12/iM. The 
efficiency e increases monotonically as the compactness gets smaller, and seems to end up with that for Schwarzschild 
black holes (e = 5.687 x 10~ 2 ). For E = 1.01// and L = 4.5/iM, the efficiency, which is e ~ 2.3 x 10 -4 , little depends 
on R/M up to the case of a black hole. 



B. Polytropic star 

From the analysis of the energy spectrum for a uniform density star, we find that the spectrum strongly depends 
on the shape of effective potential. Therefore, it may be important to see the dependence of the equation of state on 
the energy spectrum. Then, in this section, we study a polytropic star. Just for simplicity, we use the Newtonian 
polytropic equation of state as 

P = Kp l+1 ' n , (3.2) 

where n is a polytropic index and K is a constant value. 

First, we consider the first type stellar model with R > 3M. We set the polytropic index n = 1 and the constant 
K = 100 km 2 , and choose the central energy density p c = 3.0 x 10 15 g/cm 3 . This choice gives a stellar model with 
mass M = 1.267-M© and radius R = 8.862 km. The compact parameter is R/M = 4.739, and then the surface locates 
outside the peak of the Regge- Wheeler potential. We show the energy spectra of gravitational waves from a scattering 
test particle in Fig. [n]. We find the same feature as that for the uniform density star with R > 3M, i.e., there is only 
a single global peak which frequency corresponds to the orbital one of the particle at the turning point. 

Next, we discuss the second type stellar model, i.e., an ultracompact star. To construct an ultracompact star with 
a polytropic equation of state, we set the polytropic index n = 0.5, the constant K — 100 km 4 , and choose the central 
energy density p c — 388.097 x 10 15 g/cm 3 , for which the compactness parameter R/M turns out to be the smallest. 
For such a choice, the radius, mass, and compactness parameters of the star are R = 1.136 km, M = 0.296M Q , and 
R = 2.597Af, respectively. Although this seems very implausible for a realistic star, we will use this solution in our 
analysis to study the dependence of the equation of state on the energy spectrum. In Fig. [ll], we show the energy 
spectra for two trajectories of a test particle. One is the case that a particle goes into the inside of the potential 
barrier, while the other is the case that a particle passes through only outside the potential barrier. We find the same 
feature as that of the uniform density star with R < 3M . The difference is found in the frequencies of quasinormal 
modes because of the different shapes of the potential. Then, we find many small and sharp peaks in the spectrum at 
the positions different from those for a uniform density star. Consequently, we may be able to distinguish the inside 
structure of an ultracompact star by observing the energy spectrum of emitted gravitational waves. 



IV. CONCLUDING REMARKS 



We have studied the axial modes of gravitational waves from a test particle scattered by a spherically symmetric 
relativistic star. We have considered both a uniform density star and a polytropic star. We find that the energy 
spectrum depends mainly on the shape of the effective potential for gravitational waves, i.e., whether its minimum 
exists or not. 

When there is no minimum (e.g., a regular type of neutron star with a compactness parameter R/M > 3), the 
energy spectrum shows only a single global peak, which corresponds to the orbital frequency of a test particle at the 
turning point. Gravitational waves are mostly emitted around the turning point, where the effect of the gravitational 
field is strongest. 

While if a minimum of the effective potential exists (i.e., some exotic relativistic star with small compactness 
parameter R/M < 3), the energy spectrum shows a variety in its structure in addition to the global peak as the 
same as that found in the previous case. We find small and sharp periodic peaks in the region lower than w max (the 
maximum value of the effective potential), which correspond to an excitation of the axial quasinormal modes. The 
existence of the potential minimum allows the existence of a quasibound state of gravitational waves, i.e., quasinormal 
modes. Then those modes are excited by a scattered particle. This becomes more conspicuous when the particle gets 
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into the inside region beyond the potential barrier. We also find small periodic peaks in the frequency region beyond 
w max , which may come from an interference between two waves reflected by the potential barrier and the infinite wall 
near the center of the star. 

From the observation of the energy spectrum of the axial mode, what we can determine? Naively speaking, we 
may not find any information about the constituent of the star, because the axial mode does not couple to matter 
fluid. This is true for the first type stellar model, i.e., a regular neutron star. Emitted gravitational waves are mostly 
determined by the particle orbit. The case of black holes also belongs to this case, although there is a small difference 
in the high frequency part because of the difference of the boundary condition. However, if a star is compact enough 
to make a potential minimum (R < 3M), emitted gravitational waves show much more abundant information. It may 
determine a compactness parameter R/M. We also find quasinormal mode frequencies, which depend on the shape 
of the potential, and then could depend on the equation of state. 

In the axial mode case, however, in general, we cannot have direct information about the matter fluid, in particular 
about the equation of state. The polar mode of gravitational waves will provide us a more conspicuous dependence 
of matter property. Work on the polar modes is under way. 
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APPENDIX: NUMERICAL TECHNIQUES SOLVING THE LINEARIZED EINSTEIN EQUATIONS 



Here, we present numerical techniques solving the wave equations. We improved Borrelli's method p3[ to guarantee 
numerical accuracy. 

For the interior region, the numerical error of the wave function accumulates and becomes maximal at the surface, 
because the pressure vanishes there. In order to overcome this difficulty, we solve the wave equation both from the 
center and from the surface, and connect two wave functions at some inner point, e.g., r = R/2. For the exterior 
solution, we integrate the wave equation from the surface to infinity, by which we can reduce the numerical error. In 
the present model, the test particle mainly emits gravitational waves at the turning point, where we have to solve the 
wave equation very carefully because an apparent divergence appears. We shall discuss these three methods in order. 



1. Interior region of the star 



Here we show how to solve the interior wave equation (2.8). Transforming the wave function A"/™^ to Zi u = rX, 
we find the perturbation equation for Zi u as 



dr 2 



dZ, 



dr 



luj i e ^~ v 



CO 2 - ^ \ Z, ., = 0. 



(Al) 



From the regularity condition at the center, we obtain the asymptotic behavior at r ~ as 



Ziu = rn u 



J+2 



1 



2(21 + 3) 



4tt(I + 2) [-(2l-l)p c -P c ) - 



■0(t 



(A2) 



where p c , P c , and v c are density, pressure, and metric functions evaluated at the center, respectively. In order to 
reduce the numerical error at the surface, we adopt the following procedure. 

First, we construct the wave function zj® by integrating Eq. (|Al| ) from the center with the asymptotic behavior, 
Eq. flA2| ), to a middle point (r = R/2). 

Then, we find two independent wave functions zj^,zj® by integrating Eq. (Al) from the surface to the mid- 
dle point. The boundary conditions of those wave functions at the surface are given as (Z,9j , dZPy j dr) — (1,0), 

(2) (2) 

(Z\J ,dZ\J /dr) — (0, 1). By those two independent solutions, we construct the wave function Zi u (r) as 



ZUr)=a lu> z\ 1 J{r)+b lu zH\r), 



(A3) 



where the coefficients ai u and bi u are fixed by junction conditions at r = R/2. Since the wave functions constructed 
above must be continuous and smooth at r = R/2, we find the junction conditions 



R 



AD 



R 



A2) 



R 



dZ, 



(o) 



dr 



R 



dZ, 



(i) 



dr 



R i , 

77 + blu 



dZ 



(2) 



dr 



R 



(A4) 
(A5) 



From these two conditions, we find a;^ and 6;^, and then the interior wave function Z\ u 
We finally find values of X^^ and dx\™^ / dr* at the surface r = R as 



X^\R* 



Z tu (R), 



dX 



(int) 



dr* 



<"'' = I>-TT 



^ <«)->-<«)) 



(A6) 
(A7) 



where R* is the tortoise coordinate of the surface R. 
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2. Exterior region of the star 



Next we construct the exterior wave function. In order to use the Green's function method, we first consider the 
homogeneous wave function xj^ , which satisfies the wave equation 



,2 Y (° xt )(°) 



+ (u 2 - ^ (cxt) ) X^ m = 0, (A8) 



dr* 2 

and the boundary condition at the surface, which we have already evaluated in the previous subsection, that is, 

xt? (0) (R*) = xt ) (R*), (A9) 

jv-(ext)(0) jj(mt) 

— £s— (ir) = (ir) • (A10) 

Then, we construct a particular solution X^^ 1 ^ with a source term, which equation is described as 

■2 Y-(ext)(l) 



d 2 X 



^— + L 2 - X<£™ = S» (All) 



dr* 



where comes from the motion of a test particle. At the surface, the solution Xu??/ satisfies the following 

condition: 

^ K1) OR*)=0, (A12) 

,„(cxt)(l) 

— — (R*) = 0. (A13) 
dr* 

For the homogeneous equation (A8), we have two independent solutions u[™^\r*) and w[^(r*), which correspond 

to outgoing and incoming waves at infinity, respectively. Then, we describe the wave function X^^ ^ by these two 
independent solutions as 

*£? ( V) = <Wi£?(f) + A™-uS£(r*), (A14) 



where ai rnu and ^ are complex constants, which are fixed by Eqs. (A9), (A10). 

The Green's function G(r*,s*) for Eq. ( All ) is also constructed by u\™^ and as 

G( r *,s*) = ^ -u(° ut) (r*)u( in )(s*) + M< in) (r*)u(° ut )(s*) 8(r* - s*) (r* > s*), (A15) 

where W is the Wronskian and 9{x) is the Heaviside function. The wave function X/f^* is obtained by integration 

as 

/>oo 

4£ )( V)=/ G{r*,s*)S\^{s*)ds*. (A16) 



Then, the solution for the exterior wave equation (2.21) with the same boundary conditions as Eqs. (AS), (A1C) is 
given as 

r -( ext )f^*\ _ „,. _i_ o, „,, /■_*\ i / D *\o( ext ) 



^(*-*)=a jTOW « (out) (r*)+Am W «(in)(»-*)+ / G(r*,s*)SlZ'(s*)ds*. (A17) 

The asymptotic behavior of Xj^j(r*) as r — > oo is given by 

^'(r») -» (a imLJ + /w)e iwr * + (A ma , + ^K^* (r* -> oo), (A18) 
where p; mw and cr; mw are defined by 
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Plmu, w 



1 f°° 

Wr ^l^)Sl^(s*)ds*, (A19) 



a* 

oo 

(out)/ „(ext) / *n , * 



= ^ u^(0^ ; (0^. (A20) 
Since there is no incoming wave, we find the boundary condition such that 

Plmu, + Vlrrw = 0. (A21) 

Then, the amplitude Ai mw of gravitational waves at infinity is given by 

Mmuj = Otlmu + pirnuj- (A22) 

In order to extract the amplitude Ai mtJ from our solution, we use the relation 



Plrnuj Plrau ~\~ Qlmuj&lmuj — 



R* 



? (ext) , 



IS 



1 r 00 

_L_ I X^\snS^(s*W. (A23) 



Irauj V / Imuj 

R* 



With Eqs. (A21) and (A23), we finally obtain the amplitude as 



A lm ^ = a lmu + PlmuJ = - — 4? I X^ { °\s*)St2(s*)ds*, (A24) 



1 1 

Plmu, W J R « 

where Pimuj is already obtained up to an unknown constant rji u , which is canceled out with that of ^ in Eq. 



( ]A2j ) 



3. Numerical method to integrate near the turning point of a test particle orbit 

When a test particle is scattered by a star, the orbit includes a turning point (i? m in or ^min by the tortoise 



coordinate), where the radial velocity dR/dr vanishes. The source term Sj^f in Eq. (A24) is given by Eq. (2.24) 



where Qj m (u;, r) and Di m (uj,r), defined by Eqs. (2.25) and ( 2.26| ), for the present orbit (O = n/2) are described as 



/ 2 [A f 00 dR 

Q lm (u, r) = -Q m y z(I + 1)r(1 _ 2M/r) dr^6(r - R(r)) 

L_dIU(*\ i(u,T-m*) (A25) 

D lm (^ r) = ;™° lm u P drS(r - fl(r))^% (£) (A26) 



where 



r , n (m+H)/a 4 / (gj + 1) (j ~ W (V , r . 

Cim " ( X) V 47r + H)i' ( J 

which is a normalization constant of a spherical harmonics Y"; m , and P; m (cos(9) is the associated Legendre function. 



Integration over the proper time r gives the term of (dR/dr) 1 in Eq. (A26), which diverges at the turning point, 
and then the integration 

X& m P)sB{8')da' = 4= I X^°\s*)S^(s*)ds* (A28) 



w J R . lmu lmw W j j ■ 

diverges apparently. However, this is not a physical divergence, but rather comes from the change of integration 
variables from the proper time r to the radius r. This transformation becomes singular at the turning point because 
dR/dr = . In order to avoid this difficulty, we shall use a time coordinate T in the integration. 
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Since the trajectory before and after the turning point is symmetric, we choose the ini tial p osition of the particle 
T = 0, $ = at the turning point and integrate from there to infinity. The integral, Eq. ( A28|) , is then written as 



1 

W 



* i^(ext)(0) Qf(cxt) / * 



ds*X 



(**) 



°lm ~ jj. 



1(1 + 1) 



1 - 



1 

= wj 

+R (l - 

2M 
R 



(IT 
2M 



XLT } (T) S%(T) + S^(T) 



(o). 



1 ^v( cxt )(°) , s 



R J dR/dT dT 



(A29) 



2M 



S; 



(1) 



Z(Z + l)(Z-l)(Z + 2) R 
2M 



2M 
R 



D, 



slm 



8tt / 2 
'~E y/(/ + l)(Z-l)(Z + 2) 

2 fit 



1 



= -2iC, 



lm 1 



7 



l^dPlm (IT 



D, 



/ Z(Z + 1) R 3 (l - 2M/R) de V2 
AimC'im fiL 2 dP im /tt 



- ) sin^T - m$), 



(77) cos(wT - m$), 



v /2Z(r+l)(r= l)(Z + 2) i? 4 d6 
where X^ )(0) (f) = X^ m {R{f)) and dX^ m /df is its time derivative. However, the term in Eq. (|A29|) 



1 



(ext)(0) 
Imu) 



dX 



(ext)(0) 



(A30) 
(A31) 
(A32) 
(A33) 

(A34) 



is still apparently divergent because dR/dT vanishes at the turning point 

r (cxt)(0) 

we find 



To avoid this difficulty, we expand X^ ){0) , dX™ v, /dR near the turning point. From the equation of motion, 



dR 
dT 



x 



(R — -Rmin) 1 ^ 2 + 0(R — R m in), 



which is integrated as 

R-R min = a f 2 + O(f 4 ), 
where ao is an integration constant. Then, we find the wave function near the turning point as 



(A35) 
(A36) 



X 



(ext)(0) 
Ivauj 



.(ext)(0) 
" hnu) 

y (ext)(0)/A x 



dX 



(oxt)(0) 



«0 



dX 



^min 

(oxt)(0) 



Imu) 



dx, 



(ext)(0) 
Imtu 



dx 



(ext)(0) 
Imtu 



dR 



dr 



Rn 



1 - 2M/R min 

, 2y (cxt)(0) 



dr 2 



dr* 



Ru 



(A37) 



(R — R„ 



1 



dX 



(ext)(0) 



1 - 2M/R m 
a 



dr* \ mil 7 



2M dAT 



(ext)(0) 
Imtu 



(1 - 2Af/i? n 



dr* 



J2 Y ( ext )(°) 

dr* 2 



(A38) 



(ext)(0) 



dX]^. , / dr* , and d 2 Xl^ >{v> /dr* 2 at the turning point -R* nin are foun d easily by integration of Eq. 



, 2 v (ext)(0) 



whe re X t 

(A8). Since those functions are regular at the turning point, we can evaluate Eq. (A29) without any difficulty. 
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FIG. 1. Effective potentials of a uniform density star for the I = 2 mode. Solid, dotted, and dashed lines correspond to 
stellar models with radius R = 2.26M, 2.3M, and 5.0M, respectively. For an ultracompact star (R < 3M), the potential has a 
minimum, which traps gravitational waves emitted by a test particle. 
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FIG. 2. The energy spectrum (1 = 2 mode) of gravitational waves from a test particle with E — 1.01// scattered by a uniform 
density star with R = 5.0M. Solid, dotted, and dashed lines show the cases of the angular momentum L — 4.5/iM, 5.0/iM, 
and 5.5/tM, respectively. Longitudinal lines represent the frequency of a test particle at the turning point. Peaks correspond 
to orbital frequencies at the turning point of a test particle. 
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FIG. 3. The energy spectrum (1 = 2 mode) of gravitational waves for a uniform density star with R = 5.0M and for a 
Schwarzschild black hole. We choose the energy and angular momentum of the test particle as (E, L) = (1.01/x, 4.5/iM). Solid, 
dotted, and dashed lines correspond to the cases of the uniform density star with R = 5.0M, 4.0M, and of a Schwarzschild 
black hole, respectively. The peak in the spectrum coincides each other because of the same trajectory of the test particle. 
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FIG. 4. Waveforms observed in the direction of 6 = 7r/2, <f> = (I = 2 mode) for a uniform density star with R = 5.0M 
(solid line) and for a Schwarzschild black hole (dotted line) . We set the energy and angular momentum of the test particle as 
(E, L) = (l.Ol/i, 4.5/iM). In this case, the + mode cancels out in the waveform. 
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FIG. 5. The energy spectrum (1 — 2 mode) of gravitational waves for the case of a uniform density star with R — 2.26M. 
Solid and dotted lines correspond to (E,L) = (2.38/U, 12.0/iM) and (1.01/n, 4.5/xM), respectively, (b) is an enlargement of (a) 
to show the resonant peaks in detail. Longitudinal lines represent the frequency of a test particle at the turning point. We find 
some periodic peaks, which never appear in Fig. 2. Those small periodic peaks are classified into two types: one for lu < w max 
and the other for u > tJ max , where o; m ax corresponds to the maximum energy of the effective potential in Fig. 1. 



-0.03 




(Re) Mm 



FIG. 6. The energy spectrum (1 — 2 mode) of gravitational waves and axial "quasinormal" modes for the case of a uniform 
density star with R = 2.26M. Solid and dotted lines correspond to (E,L) = (2.38/1, 12.0/tiM) and (1.01/U, 4.5/xM), respectively, 
which are the same those in Fig. 5. A longitude dash line represents the frequency of the maximum w ma x of the effective 
potential. Solid (u < w ma x) and open (to > u max ) circles denote axial "quasinormal" modes of the uniform density star with 
R = 2.26M. We find that resonant periodic peaks for ui < w max agree quite well with those quasinormal modes, which are 
excited by a scattered test particle. However, small periodic peaks for to > w m ax seem not coincide exactly with "quasinormal" 
modes. 
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FIG. 7. The waveform observed in the direction of 9 — n/2, <f> = (I — 2 mode) for a uniform density star with R — 2.26M 
(solid line) and for a Schwarzschild black hole (dotted line). We set the energy and angular momentum of the test particle as 
(E, L) = (2.38//, 12.0/iM). In this case, the + mode cancels out in the waveform. 
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FIG. 8. The energy spectrum (1 = 2 mode) of gravitational waves for uniform density stars with two different compactness 
parameters, R = 2.26M (solid line) and 2.3M (dotted line). We set the energy and angular momentum of the test particle as 
(a) (E,L) — (2.38/i, 12.0/iM) and (b) (E,L) = (1.01/n, 4.5/iM). Although global peaks coincide with each other because of the 
same trajectory of the particle, resonant periodic peaks are different because of the different compactness. This may show us 
how to distinguish the compactness of a star by use of observational data. 
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(a) (b) 
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FIG. 9. The energy spectra (1 = 2 mode) of gravitational waves for a uniform density star with R = 2.26M (solid 
line) and for a Schwarzschild black hole (dotted line). We set the energy and angular momentum of the test particle as (a) 
(E,L) = (2.38/i, 12.0/j.M) and (b) (E,L) = (1.01/U, 4.5/nM). This difference may enable us to distinguish an ultracompact star 
from a black hole. 




FIG. 10. The energy spectrum (1 = 2 mode) of gravitational waves for the case of the polytropic star with a large 
compactness parameter. We set the polytropic index n = 1, the center of the density p c = 3.0 x 10 15 g / cm 3 , and the energy 
of the test particle as E = 1.01/4. The solid, dotted, and dashed lines correspond to the angular momentum of the test particle 
L — 4.5/iM, 5.0/xM, and 5.5/ttM, respectively. Longitudinal lines represent the frequency of a test particle at the turning point. 
This figure shows the same feature as Fig. 2. 
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FIG. 11. The energy spectrum (1 = 2 mode) of the gravitational waves for the case of the polytropic star with a small com- 
pactness parameter. We set the polytropic index n = 0.5 and the center of the density p c = 388.097 x 10 15 g / cm 3 , which may be 
implausible but is necessary to obtain an ultracompact star. The solid and dotted lines correspond to (E, L) — (2.38/i, 12.0/iM) 
and (1.01//, 4.5/iM), respectively. This figure shows the same features as Fig. 5. 
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TABLE I. Axial "quasinormal" modes for a uniform density star with R = 2.26M. We use the continued fraction expansion 
method, which was first used to calculate black hole quasinormal modes by Leaver [24] , and w modes of a polytropic star by 
Leins, Nollert, and Soffel [16] . Our results agree quite well with those of Kokkotas, who adopted a different method [20]. Our 
imaginary part has opposite sign to Kokkotas's result, because our definition of the Fourier transformation is opposite to his. 
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